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Abstract. In this article we prove L p estimates for resolvents of La- 
place-Beltrami operators on compact Riemannian manifolds, generaliz- 
ing results of [12] in the Euclidean case and [17] for the torus. We follow 
|18j and construct Hadamard's parametrix, then use classical bounded- 
ness results on integral operators with oscillatory kernels related to the 
Carleson and Sjolin condition. Our initial motivation was to obtain IP 
Carleman estimates with limiting Carleman weights generalizing those of 
Jerison and Kenig [TT]; we illustrate the pertinence of L p resolvent esti- 
mates by showing the relation with Carleman estimates. Such estimates 
are useful in the construction of complex geometrical optics solutions to 
the Schrodinger equation with unbounded potentials, an essential device 
for solving anisotropic inverse problems [4]. 



This article aims at proving LP estimates on the resolvent of Laplace- 
Beltrami operators on compact Riemannian manifolds in the spirit of those 
obtained by Kenig, Ruiz and Sogge |12j on the flat Euclidean space and of 
Shen [17J on the flat torus. This work grew out of our recent interest in 
LP Carleman estimates with limiting Carleman weights generalizing Jerison 
and Kenig's estimates [TT] to variable coefficients and used in the context 
of anisotropic inverse problems [HE]- It turned out that in [S], we were 
able to prove LP Carleman estimates in a direct way, without referring to 
resolvent estimates. Meanwhile, in the process of obtaining those Carleman 
estimates, we were led to derive resolvent estimates. Since these resolvent 
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estimates might be helpful in different contexts, we decided to make them 
available. 

Let (M, g) be a compact Riemannian manifold of dimension n > 3 (with- 
out boundary). In local coordinates, the Laplace-Beltrami operator takes 
the form 

^ g y/det g dxj i^^ 99 dxt, 

where (<r )i<y,fc< n denotes the inverse matrix of (gjk)i<j,k<n if 9 = 9jk dx- 7 ® 
dx k . As is classical, we are using Einstein's summation convention: repeated 
upper and lower indices are implicitly summed over all their possible values. 
We will be using the notation 

\0\ 9 = \l 9*016* 

to denote the norm of a vector and d g to denote the geodesic distance in 
(M,g). We denote er(— A g ) the set of eigenvalues of the Laplace-Beltrami 
operator on (M,g), and denote by 

R(z) = (-A g + zy\ z€C\a(-A g ) 

the resolvent operator. We are interested in proving estimates of the form 

\W\\l\m) < C\\{A g - z)u\\ L a {M) 

for the couple of dual exponents (a = 2n/(n + 2),b = 2n/(n — 2)). Note 
that by taking u = 1, it is clear that such estimates can only hold away from 
z = 0, and similarly by taking u to be any eigenfunction of — A g , away from 
the negative of any eigenvalue. This is in contrast to the non-compact flat 
case where Kenig, Ruiz and Sogge |12j were able to obtain such estimates for 
all values of the complex parameter. In fact, we will restrict our investigation 
to the following set of values 

E 5 = {z e C\R_ : Rev^> S} 

= {zeC\R_: (Imz) 2 > A5 2 (5 2 - Rez)}. 

This set is the exterior of a parabola and can be visualized in figure [TJ The 
corresponding results were obtained by Shen [T7] on the flat torus. 

It seems that such estimates cannot be obtained as a direct consequence of 
the spectral cluster estimates obtained by Sogge [18]. Indeed, if Xm denotes 
the cluster of spectral projectors related to the eigenvalues of the Laplacian 
the square root of the negative of which lie in the interval [m, m + 1) then 
Sogge's estimates read 

\\ Xm u\\ 2n < C{l+m)^\\u\\ L 2 {M) 

Ln-2 (M) V ; 

\\XmU\\L 2 (M) < C(l + m )' 2 \\ u W T ^ 7 r„,- 

This would reduce the resolvent estimate to an L 2 estimate of the form 

IIXm^||L 2 (Af) < A ( m , z )\\Xm( A 9 ~ z HIl 2 (M) 
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but summing the corresponding series ^™ =0 (l + m)A(m, z) to patch the es- 
timates requires a decay of A(m, z) which cannot be achieved. Nevertheless, 
we will use Sogge's estimates in the process of reducing Carleman estimates 
with limiting Carleman weights to resolvent estimates in the last section of 
this paper. 




Figure 1. Allowed values of the spectral parameter z 



Theorem 1.1. Let (M,g) be a compact Riemannian manifold (without 
boundary) of dimension n > 3, and let 5 £ (0,1) be a positive number. 
There exists a constant C > such that for all u £ C°°(M) and all z £ Eg, 
the following resolvent estimate holds 

(1.1) ||U|| 2n < C||(A ~ Z)U\\ 2n 

The proof of the theorem relies on the construction of the Hadamard para- 
metrix. We provide a fairly complete study of the LP boundedness properties 
of Hadamard's parametrix and the corresponding remainder term. Some of 
the estimates might not be used in the proof of Theorem |1.1| but may be of 
independent interest in other contexts. 

Remark 1.2. The question of whether the estimates hold for a larger set 
of values of the spectral parameter z, for instance 

Sj = {z£C \ R : | Im z\ > 5 when Re z < and \z\ > 5 when Re z > 0} 
remains open. 

As is common practice, we will write A < B as shorthand for A < CB 
where C is a constant which depends on known fixed parameters (the di- 
mension n, the indices p, q of the Lebesgue classes involved, the Riemannian 
manifold (M,g), etc.), but whose value is not crucial to the problem at 
hand. Similarly we use the notation A ~ B when A < B and A > B. We 
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consider the principal branch of the square root on C \ R_ which we denote 
by yfz. We also choose the convention that arg z E (— it, tt) when z £ C\R_. 

In the second section of this note, we review the different bounds known on 
oscillatory integral operators, in particular those obtained by the Carleson- 
Sjolin theory. In the third section, we construct the Hadamard parametrix 
following [9j. In the fourth section, we prove bounds on the Hadamard 
parametrix following the techniques of |12l [18] and derive the corresponding 
resolvent estimates. Finally, in the last section, we illustrate the pertinence 
of resolvent estimates by showing how they imply Carleman estimates with 
limiting Carleman weights in the spirit of Jerison and Kenig's estimates |llj . 
Note that these Carleman estimates can be obtained by other means 

Acknowledgements. Carlos Kenig was partially supported by NSF grant 
DMS-0968472. Mikko Salo is supported in part by the Academy of Finland. 
David Dos Santos Ferreira and Mikko Salo would like to acknowledge the 
hospitality of the University of Chicago. 

2. Oscillatory integral operators 

This section is a short review of the classical boundedness results of oper- 
ators with oscillatory kernels which will be needed in this note. We consider 
operators of the form 

T x u(x)= [ e a ^a(x,y)u(y)dy 



where a belongs to Cq°(C7 x V) with U, V open sets in R™ and <p is a smooth 
function on U x V. The product T*U x T*V of cotangent bundles of U 
and V endowed with the symplectic form a = d£ A dx + dv A dy is a sym- 
plectic manifold of dimension An. A fundamental object in the study of the 
operators under our scope is the Lagrangian submanifold of T*U x T*V 

= {(x,d x ip(x,y),y,d y <p(x,y)) : (x,y) G U x V}. 

The first result relates to L 2 boundedness of such oscillatory integral oper- 
ators. For that purpose, we consider the two projections 

C v C T*U x T*V 



T*U T*V. 
If 7T{7 is a local diffeomorphism, i.e. 

(2.1) det— -£-^0 onUxV 

oxoy 



then so is 7iy and C v is locally the graph of the canonical transformation 
otl 

to the case where C v is locally a canonical graph. 



q = Tr^oTry. The mother of I? estimates [5] |2H Theorem 2.1.1] corresponds 
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Theorem 2.1. Suppose that the non- degeneracy assumption (2.1) is satis- 
fied, then for all compact K C V there exists a constant C > such that for 
all A > 1 and all u £ C^{K) 

\\T\u\\ L 2 < C\- n/2 \\u\\ L 2. 

This non-degeneracy assumption is too stringent for the applications in 
this note (we stated this result only for the sake of completeness), in fact, we 
will use a weaker result and make the following non-degeneracy assumption 
Theorem 25.3.8] 

(2.2) corankcr^ = 2 
on the lifted symplectic form 

Note that we could have used either the first or the second projection ny to 
lift the symplectic form; the condition is explicitly given bj£] 

(2.3) corank — = 1 on U X V. 

oxoy 

The L 2 boundedness result that we will need reads as follows. 



Theorem 2.2. Suppose that the non- degeneracy assumption (2.3) is satis- 
fied, then for all compact K C V there exists a constant C > such that for 
all A > 1 and all u G C^(K) 

\\T x u\\ L 2 < C\- {n - l ^ 2 \\u\\ L 2. 

By interpolation with the trivial L°° — L 1 bound, one gets the following 
LP - LP' result. 



Corollary 2.3. Under the assumptions of Theorem 2.2, we have 
\\T x u\\ LP , <C p \-( n -V/v'\\u\\ LP 

where 1 < p < 2. 

The range of admissible exponents in this theorem is represented by the 
green line in figure [2} 

The second type of results we need are those related to the n x n Carleson- 
Sjolin hypothesis pQ. Under this curvature assumption, the former L p — L q 
boundedness results may be improved by a factor A -1 / p . We refer to |22|, 
Chapter IX] and [2H Chapter 2] for a more detailed exposition of such 
results. Let us start by recalling this assumption: suppose that (2.2) holds, 
then 

EV = {d x <p(xo,y):yeV}cT* U 

Xl = {d y <p(x,yo):xeU}cTZ V 
are hypersurfaces in the respective cotangent spaces. 

^This is also equivalent to the fact that dnu : T{C V ) —¥ T(T*U) has rank In — 1. 
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Figure 2. Admissible exponents in Theorem 2.4 and Corol 



lary 2.3 



Carleson-Sjolin Assumption. Both hypersurfaces and T>y Q have ev- 
erywhere non-vanishing Gaussian curvature for all xq G U and yo G V . 

Under the nxn Carleson-Sjolin hypothesis, we have the following bound- 
edness result [21, Corollary 2.2.3]. 



Theorem 2.4. Suppose that the non- degeneracy assumption (2.3) is satis 



fied and suppose that the nxn Carleson-Sjolin condition holds. Then for 
any compact set K C V there exists a constant C > such that for all A > 1 
and allue C^(K) 

(2-4) ||Tau|U <( r») < C\- n / q \\u\\ LP{nn) 

with q = andl<p<2. 

The range of admissible exponents in this theorem is represented by the 
red line in figure [2} 



Let us finish by the following important remarks. 

Remark 2.5. If the amplitude a depends on parameters, the estimates in 
the theorems are uniform with respect to the parameters as long as the 
amplitude is supported in a fixed set XJ x V and is uniformly bounded as 
well as all its derivatives. Similarly if the phase if depends on parameters, 
the estimates remain uniform as long as all its derivatives are uniformly 
bounded and either 

- there is a uniform lower bound on a minor of the determinant of the 



mixed Hessian in (2.3) in the case of Theorem 2.2 and Corollary 2.3 
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or there is a uniform lower bound on the Gaussian curvature of the 



hypersurfaces in the case of Theorem 2.4 



Remark 2.6. Although Theorems |2.3| and |2.4| are stated for oscillatory 
integral operators T\ with amplitudes a G C^(U x V) with constants that 
depend on the compact set K = suppa, they are still valid with constant 
depending on R, A where we assume that 

{x — y : (x, y) G supp a} C B(0, R) 

that a and all its derivatives are bounded on U x V with bounds only 
depending on A and with phase (p satisfying non-degeneracy or curvature 



assumptions of the type discussed in Remark 2.5, uniformly on U xV, with 
constant depending only on 1/A. 

By the compact support assumption on the amplitude a in U x V , we have 
a natural extension of T\ to an operator acting on functions in R™, which 
gives a function on R™. We proceed as in [221 pages 392-393]. Firstly, one 
observes that it suffices to prove an estimate of the form 

( 2 - 5 ) \\T\ u \\li(B(x ,1)) < C^~ S \\u\\LP(B(x ,r)) 

with a constant C uniform with respect to xq G R n and some radius r 



depending only on R. In fact, if (2.5) holds, raising both sides of (2.5) to 



the t/-th power q and integrating with respect to xq G R ra gives 

|B(0,1)| \\T x u\\ q Lg < C q \- sq [ ( [ \u(x)\ p dx) P dx 

J \J\x-x \<r J 

and by Minkowski integral inequality (q > p) 



|fl(0,l)|||r A ti||£, <C q \~ sq (^J \u(x)\ p \B(Q,r)\« dx 



so that 



Secondly, the estimate ( |2.5| ) follow from the LP ' — L q Theorems 2.3 and 2.4 
discussed in Remark |2.5| because the operator T\ does not move too much 
the support: let x G Cg°(B(0, 2)), resp. V G C$°(B(Q, 2r) be cutoff functions 
which equal one on B(0, 1), resp. B(0,r). Then one has x( x ~ x o) a { x ^y) = 
x(x — Xo)a(x, y)ijj(y — yo) for r large depending on R and 

11^11x3(^(20,1)) < ||x(- - xo)T\u\\ L i = \\x(- ~ x )T x ip(- - x )u\\ Lq . 



Finally, the estimate (2.5) follows from the LP — LP boundedness results on 



the oscillatory integral operator 

X (x)a(x + x ,y + y )iP(y)e a ^ x+X0 'y + y^u(y) dy 
which are uniform in xq by our assumptions on a, (p. 
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3. Hadamard's PARAMETRIX 

We follow [18] . and use Hadamard's parametrix (9j Section 17.4] to prove 
the resolvent estimates. We first introduce the following sequence of func- 
tions 

including a fundamental solution Fq of (— A+z) where A is the flat Laplacian 
on R™. We used | • | for the Euclidean norm (recall that | • \ g denotes the 
norm induced by the Riemannian metric). These are radial functions and 
in polar coordinates, the equation on F v reads 

/ » \ n — 1 _ „ \ Sn when v = 

-A + z)F v = -8 2 T F U d r F u + zF v ={ 

r I uF v -i when v > 

Besides, as can be seen from their Fourier transform, these functions satisfy 
(3.1) 8 r F v = ~\Fv-\, v > 0. 

3.1. Construction of the parametrix. Let xo G M and U be a neigh- 
bourhood of xq in M. Since the construction is local, we can think of U as 
a ball in R ra , and work in coordinates. If U is small enough, the distance 
function r = d g (x,y) is smooth on the product U x U with the diagonal 
removed, and furthermore one can take polar normal coordinates for x in U 
with center at y. In polar coordinates 

x = exp y (r6), 6eT y M 

the volume form takes the form 

dV g = r n - 1 J(r, 9) drAdO 

where d# denotes the canonical measure on the unit sphere S y M of the 
tangent space (T y M, g(y)). The Laplacian of the distance function r = 
d g {x,y) is given bjQ 

n-1 d r J(r,9) 
9 r J(r,6) 

and this implies for a radial function 

A g f(r) = f"(r)\dr\ 2 g + f'(r)A g r 

= f(r) + ^V(r) + dj j-f{r). 

This identity holds in the classical sense when x ^ y and in the distribution 
sense on U. If we use those computations on the function ooFq, we get 

(-A g + z)(a F ) = a 5 - 2^d r a + ^J a o^9 r F - (A g a )F 
2 We refer the reader to [B] section 4.B] for those computations. 
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-A g + z)F = a S(r) - 2{ d r a + -^ao )d r F 



and on the function a u F u , v > using (3.1) 

(-A g + z)(a v F v ) = (rd T a v + (u + r^Ha^i^-i - {A g a u )F u 
in V'(U). Thus, if we set 

AT 

F(x,y,z) = ^2a u (x,y)F u (d g (x,y), z) , 

then we have 

d r J 

+ ^2 (^drOiv + (v + r ^pJ a u ~ Ag<Xv-lj F u-i ~ (A g a N )F N 

in T>'(U), where the derivatives were taken with respect to x. Here S(r) 
stands for the pullback of the Dirac mass 5q on R by the map x h-> d g (x, y), 
this can be explicitly computed 

1 

y/det g(y) 

If we choose the coefficients a u to be solutions of the equations 

d r J 

(3.2) d r a + -^j a o = 

(3.3) rd r a v + (^u + r^p S ja u = A g a u -i, l<v<N 
then the function F satisfies: 

(3.4) (-A g + z)F = a 5(r) - (A g a N )F N . 



The equation (3.2) can be explicitly solved by taking 

_ i 

i + / N-i/4 1 / detfir(y) N * 

ao = det g(y) ' = — —, — - 

\/J[r,0) Vdet(exp*g)(expg x) , 

Note that the computations that we just made are not coordinate invariant, 
neither is the quantity det g (contrary to the volume form) so the choice of 
ceo may be understood in the following manner: first pick normal coordinates 
in x depending on y 

(x,y) -> (exp-^x),^) 

then ao depends on the ratio of the determinant of the metric expressed in 
the coordinates in y and of the determinant of the metric expressed in the 
normal coordinates (depending on y\) in x that we have just picked. Having 
given those precisions, it is clear that the function ao is smooth and satisfies 
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Setting a v = a^ x a v and /3„_i = o.^ 1 A g a u -\ reduces the equation (3.3) 

to 

rd r a u + va v = 3 u -i, 1 < v < N . 

V v ' 

=r 1 - v d r (r v a v ) 

In polar coordinates, this is explicitly solved by 

a v = r~ v \ t u ~ 1 p v -i(exp te)dt= / ^ _1 /3i/-i(expj, ri0) di, i/ > 1. 
By induction, one can solve the equations on a u by taking 

Jo 

where j xy is the distance minimizing geodesic from y to x (this is well 
defined if the neighbourhood U of xq is small enough). These functions are 
smooth on U. 

Let x ^ Co°(U x ^0 be a symmetric function which equals one near the 
diagonal, we get using our computation (|3.4l) 



(3.5) (-A g + z)\x(-,y)F(-,y,z) 

+ [x(-,y),A g ]F(-,y,z)+x(-,y)(A g a N )F N (d g (-,y),z) . 



x(y,y) x 

i y 

V det 9(y) 



=H N (x,y,z) 

Note that [x(-,y),A g ] = 2gr&d g x(-,y) is a first order differential operator 
whose coefficients are supported in the support of x an d away from the 
diagonal x = y. We are ready to consider Hadamard's parametrix 

Tft ad (z)u= / x(x,y)F(x,y,z)u(y)dV g (y) 
Jm 

which satisfies 

(3.6) (-A g + z)T Ua ,d(z)u = x(x, x)u + S(z)u 

with S(z) = 2grad 9 x(-,y) °T RaA {z) + S 2 (z) = Si(z) + S 2 {z) and 



S 2 (z)=i H N (x,y,z)u(y)dV g (y). 
Jm 

The integer N will be chosen, in the next paragraph, large enough so that 
S 2 (z) is smoothing. Given the fact that the function 

N 

F(y,x,z) = ^2a„(y,x)F u (d g (x,y)) 

u=0 

behaves in a similar fashion to F(x, y, z) (the only lack of symmetry comes 
from the coefficients a v whose only relevant property is their smoothness), 
the transpose t Tn ai d( z ) shares the same boundedness properties than Tn^z). 
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It is therefore not misleading to think of the Hadamard parametrix as a sym- 
metric operator in terms of its boundedness properties. 

3.2. Bessel functions. In this paragraph we describe the behaviour of the 
functions F v (see [121 page 338-339], jTSJ Lemma 4.3]). They can be explic- 
itly computed in terms of Bessel functions 

(3.7) F v {r,z) = cr-f+^f-^ K n/2 ^_ u (^r) 

where K n i2-\- v are Bessel potentials 



K m (w) 



-w cosh t 



cosh(mt)di, Kew > 0. 



The qualitative properties of K m are as follows 

\K m (w)\ < C m \w\~ m when \w\ < 1 and Rew > 0, 

K m (w) = a m (w)w~2e~ w when \w\ > 1 and Rem > 0, 

where the functions a m have uniform estimates of the form 



(3.8) 



d / w 

T, «m f- : 

or a \ \w\ 



a e N 



when r > 1. We sum up the well known properties of the function F v 
obtained from the explicit formula (3.7) in the following lemma. 

Lemma 3.1. There is a constant C > such that 

\F v (r,z)\ < Cr- n+2+2u when r < \z\~*. 

Furthermore, one has 

F u {r,z) = \z\ S ~i k ~ E ^ L e'^r^^^ a u (r, z) 

where the functions a u have a symbolic behaviour 







— ( r ) z) 



^ C a ,v r 



a G N 



when r > \z\ 2 . The constants are all uniform with respect to z G C, 
\z\ > c. 



In particular, if N > (n — l)/2 then by Lemma 3.1 we have 
(3-9) \F N {r,z)\<\z\- l l 2 
when \z\ > 5. We also recall that 

N 



F(x,y,z) = ^2a„(x,y)F u (dg(x,y),. 



u=0 



with smooth coefficients ct v . Thus from Lemma |3.1| we deduce bounds and 
asymptotics on the function F, which we synthesize in the following lemma. 
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Lemma 3.2. There is a constant C > such that 

(3.10) \F(x,y,z)\ < Cd g (x,yy n+2 when d g (x,y) < \z\~*. 
Furthermore, one has 

(3.11) F(x,y,z) = Izl^-se-v^sfof) dgix^y)'^ a(x,y, z) 
where the function a has a symbolic behaviour 

\d% y a(x,y,z)\ < C a d g (x,y)~ a , a G N 

when d g (x,y) > \z\~z. The constants are all uniform with respect to z € C, 
\z\ > c. 

4. Resolvent estimates 
This section is devoted to the proof of LP estimates on the resolvent. 

4.1. Parametrix estimates. Our goal in this paragraph is to prove that 
the Hadamard parametrix is bounded on certain Lebesgue spaces. Due to 
the local nature of the construction, the kernel of the parametrix is supported 
in U x U with U a small neighbourhood of a point xq £ M; using local 
coordinates, we might as well suppose that xq = and that U is a small 
ball in R n of radius e > 0. 

Theorem 4.1. The Hadamard parametrix is a bounded operator 
(4.1) \z\ s T^ d {z):U\M)^L\M) 

with a norm uniform with respect to the spectral parameter z £ C, \z\ > 1 
when 0<s<l,p<2<q 

1 1 2s 2 

4.2 + - = -, 

p q n n 

and 

(I 1 1 1\ 1 2 112 

min , > — , < < — . 

\p 22 q J 2n n+1 p q n 



Remark 4.2. The region of admissible exponents in Theorem 4.1 is the 
trapezium CC'D'D (blue quadrilateral in figure [3]) with vertices 

\2 2n 2 2n)' \2 2n 2 2n 

D= (1 1 (ji-lf \ D ,_fn 2 + 4n-ll I 



2 2n'2n(n+l)7' \ 2n(n + 1) '2 2n 

Note that 

F: /i ii i 



2 n+12 n+1 



is the middle of [CC'j. 
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± ± - ^ I 

Figure 3. Admissible exponents in Theorem 4.1 



Let -00 £ Cg°(R) be supported in [—1,1] and equal one on [—5,2]) let 
ip = ipo(-/2) — tpo E C^°(R) be supported in [— 1, — |] U [i,l] 3 we consider 
the following dyadic partition of unity 



u=0 

We denote Vi/ = "0(2 _I/ -) if i> > 1 and decompose Hadamard's parametrix 
in the following way 

00 

T Had (z) =T Q (z)+J2 T »{z) 

where 

T v (z)u(x) = / x(x,y)ipu{\z\^d g (x,y))F(x,y,z)u(y)dV g (y). 
J: ' 

We first observe that 
(4.3) 



M 



1 



4(x,y) = ^G jk (x,y)(x j - y j )(x k - y k ) 



where G(x, x) = g(x) for all x G U. Indeed, the square of the distance 
function vanishes at order two at x = y and its Hessian is twice the metric 
since 



2 



t=o 



2ie 
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In particular, in the chart U the geodesic distance is equivalent to the Eu- 
clidean distance, i.e. 

(4.4) —\x-y\<d g (x,y)<C\x-y\. 
when x,y £ U. 



Let us begin by examining the term Tq(z): using the bound (3.10) in 



Lemma 3.2, we majorize the kernel by a constant times 

I i-n+2-i 

\ x -y\ l \x-y\<C~ 1 \z\- 1 / 2 - 

Using Young's inequality when < 1/p — l/q < 2/n and Hardy-Littlewood- 
Sobolev's inequality when 1/p — 1/q = 2/n, and scaling the estimates, bear- 
ing in mind the condition (14. 2|), we get 



n — 2 , n 2s il 1 i 

(4.5) \z\ s \\T (z)\\ ciLP , Lq) < \ z \-^ + -A^ + - P - q - 1 ) = 1 

when < 1/p — 1/q < 2/n excluding the endpoints (0,2/n) and (2/n, 0) of 
the segment 1/p — 1/q = 2/n. Now we deal with the other terms: on the 
support of the kernel of T v (z) we have d g (x,y) > 2 u ~ 1 \z\~ 1 / 2 and if we use 



the asymptotic description (3.11) in Lemma |3.2[ the operator T u (z) assumes 
the form 

T u (z)u(x) = 

Im d g (x,y) 2 



n-l 1 / 

\z 4 2 / e 



If we denote 



e -Re^zd g (x,y) 



A(x,y,z) = —j-g(x, y, z)x{x, y)y / detg{y) 

d g (x,y) 2 

then we have 

(4.6) d« y A = o(d g (x,y)-^-W 

uniformly in z. 

We need to distinguish two further cases: argz G [— ir/2,ir/2] where the 
exponential is uniformly decaying and argz ^ [— tt/2,tt/2] where the expo- 
nential has an oscillating behaviour. Let us begin by the first case which is 



straightforward: thanks to (4.4) the kernel of T u (z) may be bounded by a 
constant times 

i ,11=2 — £=|^| 1 /2| :I ._ 2/ | 

\Z 2 2 2 V Q V2 1 1 1 Wl 



and Young's inequality provides the bound 

\c(LP,Li) 



z\ s \\T u {z)\\ c[Lv r q) <2-^e- c2V . 
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This (together with (4.5)) can be summed as a geometric series 



oo 



A" \\ T B a d{z)\\c(LP,Li) < \ Z \ S W T v( Z )\\jC(LP,L9) 



oo 



<Vr^<i 



so it only remains to consider the case where argz ^ [— ir/2, vr/2]. 



The case - — - > 5 + 2n • We majorize the kernel by a constant times 




n — 2 n — 1 



\ Z \ 2 Z 2 1 \x-y\<c- 1 2"\z\- 1 / 2 

and using Young's inequality we get 



z\lMz)\\ c(LP , Lq) <2^ +n (l~ 1 P )). 



As above we can sum a geometric series and obtain the boundedness of the 
M s (^Had(-20 — Tq(z)) as long as the exponent is positive, i.e. when 

1111 

p q 2 2n 

This region is in the plane (1/p, l/(/) the triangle ABB' (delimited by the 
purple lines in figure [3]). 

The other cases. We first observe that, because of the support properties 
of the kernel, the only terms T u (z) which are nonzero are those for which 
2 u \z\~ 1 / 2 is bounded, and in fact even small if one is prepared to shrink the 
size of U. We consider the following rescaled distance function 




~ \x — y 



and the following rescaled operator 





which was obtained by the 
phase is given by 




and the amplitude by 
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With those notations, we have 

T v {z)u(x) = {2- y \z\^y n f v {z)v(2- v \z\^x), v(x) = u(2 v \z\~*x). 

Because of ( |4.6[ ) and because 2 u \z\~ 1 / 2 < 1, both phase and amplitude are 
uniformly bounded as well as all their derivatives. Besides, the amplitude 
a u (-,-,z) is compactly supported in R n x R n (with possibly a very large 
support) but the set 

{x-y: (x,y) £ suppa„(-, -,z)} 

is contained in a fixed compact set of R n . 

This is an oscillatory integral operator, whose phase is (a rescaled version 
of) the Riemannian geodesic distance, and the previous observation shows 



that part of the assumptions in Remark 2.6 are fulfilled. Besides we have 



d 2 ip u ( arg z\ d 2 

sin 



dxdy \ 2 J dxdy 
v v ' 

I— |>V2/2 



9jk(x )(xi -yi){x k -y k ) + 0{e) 



this means that the non-degeneracy condition (2.3) is satisfied if e is taken 



small enough. Corollary 2.3 applies 

(4.7) \\f u (z)u\\ Lp/{M) < \z\^2-^2-^- 1 ^'\\u\\ LP{M) 
and after rescaling, we obtain 

\*\'\\T»{z)\\cW) $ (\z\ 1/2 2-n n ^7- 1 )\z\^2-^2-^- 1 W 

,.( n + l n-l \ 

(4.8) <2"y—-—). 

Again, we can sum those estimates into a geometric series provided that 

1 1 1 

> r + 



p 2 n + l 

This shows that the operator |z| s (TH a d(£) — Tq(z)) is uniformly bounded on 
the semi-open segment [AF) in figure [3j 

Since the distance function r = d g (-,y) satisfies the eikonal equation 

|dr| s = 1 

its differential parametrizes a piece of the cosphere bundle, thus 

(4.9) S yo = {d x d g (x, y ),xGU}c S* Q M 

has everywhere non- vanishing Gaussian curvature. This is also true when 
x and y are interchanged because of the symmetry of the function. This 
remains true after rescaling and it is not hard to see that there is a uni- 
form lower bound on the Gaussian curvature when e is small enough. This 
means that the Carleson-Sjolin condition is satisfied (in its uniform variant 



mentioned in Remark 2.6) and Corollary 2.4 applies 



(4.10) \\T u (z)u\\ Lq{M) < \z\^2-^2-^\\u\\ LP{M) 
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if q = (n + l)p' /(n — 1) and 1 < p < 2. After rescaling we finally obtain 

\z\ s \\T„(z)\\ c{LP;Lq) < (Izl^JTT^--- 1 ) 1*1^2-^2-^ 
(4.11) < 2 -Hf-— ). 

Those estimates can be summed into a geometric series when 

1 1 1 

p > 2 + 2n" 

This shows that the Hadamard parametrix is bounded on the red segment 
[AC) in figure § 

If we sum up the results obtained so far: we have shown the boundedness 
of \z\ s (Tn^diz) — Tq(z)) on the segments [AS) (by Young's inequality), [AD) 



(by the Carleson-Sjolin theory), [AF) (by the non-degeneracy corollary 2.3) 
in figure[3| By interpolation and duality, this shows that |^| s (TH a d(-2) — Tq(z)) 
is bounded on the pentagon ABDDS' in figure [3| This ends the case 
argz ^ [— 7r/2,7r/2] and hence the boundedness of Hadamard's parametrix 



(4.1). This completes the proof of Theorem 4.1 



Remark 4.3. In fact, the computation made to obtain the bound on the 
segment [A, F) is redundant since the bound can be obtained by interpola- 
tion and duality from the bound on the segment [A,D). 

4.2. Additional bounds. If one is prepared to loose some powers of the 
spectral paramater \z\, one can improve the range of allowed exponents. 



Indeed, the estimate (4.11) can still be summed on 2 U < |z| 1//2 when | < 



1 + - 1 - 

2 ' 2n 



\z\ s \\Tn a d(z)\\c(LP,Li) < ^2 W \\T„(z)\\c(lp,li) 



u=0 



E„ ( n + l n\ ( n+1 n \ 

2 U \— — p) < wn — 
rs_/ I I 

v>0:2 v \z\-V 2 <l 



Note that when q = ^j- p' we have 



(n + 1) = —-(n-1) 



P Q / P 



and therefore 



n-1 I 

\z\ 4 \p q ' 2 . 



The values of (1/p, 1/?) obtained in this estimate correspond to the semi- 



open segment (DE] on the figure pi Similarly, the estimate (4.8) can still 
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Figure 4. Admissible exponents in Theorem 4.4 



be summed when - < ^ + 

oo 

\z\ s \\T Uad (z)\\ c{LPLp/) <^\z\ s \\T,(z) 



\C(LP,LP') 



u=0 



EJ n+l n-X \ I n + l 

rs_/ I I 

i/>0:2 ! '|«|- 1 / 2 <1 



Note that we have 



(n+l) 



1 1 



2(n + l) , 



p pi J p 
and therefore we retrieve the same bound as above 



n+l n-1 \ n — 1 (1 1) 1 

. 2p' 4 / = l^l 4 Vj 9/ 2_ 



The values of (1/p, 1/q) obtained in this estimate correspond to the semi- 
open segment (FG] on the figure [61 



By duality and interpolation, we get the following theorem. 
Theorem 4.4. The Hadamard parametrix is a bounded operator 
(4.12) T Ha dW : L p (M) -> L^M) 



when 1 < p < 2 < q and 

1 1 2 

< 



n-1 1 1 n+11 

■ < - < — ■ 

p q n + l n+l p' q n — 1 p' 
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with a norm bounded by 

n-l I 1 l\ 1 

(4.13) IIW*) \\ c{LPiLq) <C\z\ — y~ P *) 2. 

The extended set of admissible exponents is the pentagon DEGE' D' , rep- 
resented by the union of the green and red sets in figure [6j These estimates 
(or rather the endpoint p = 2, q = 2(n + l)/(n— 1)) were proved by Sogge 
in [H] 

||2Had(*)w|| 2(n+i) <\z\ 4( " +1) \\u\\ L 2 (M) . 

L »-l (M) ^ ; 




1 

n-l 



1 

n+1. 



1 

n+1 n+1 

Figure 5. Decay of the parametrix: ||THad(^) l *||/:(i>,i«) = ®i\ z \ 



_2 

n+1 



The combination of Theorems 4.1 and 4.4 gives the following bound on the 
parametrix 

(4.14) \\T^{z)u\\ c{LP>Lq) < C\z\-° 

in the polygon CDEGE'D'C (except for the segments [CD], [DD'\ and 
[D'C]), where the order a is a piecewise linear function of d = 1/p — 1/q 

n- 1 , 1 , „ , 2 



(4.15) 



a 



71 



1 , 1 

^ d+ 2 
d+ 1 



when < <i < 
2 



when 



n + 1 

< d < 1 



2 n + 1 

The graph of this function is represented in figure [5} 

We finish by observing that as a function of the y variable, the kernel of 
the 7H a d(-z) — Tq(z) is L' r -integrable when r < 2n/(n — 1) since 



\(Tn 



ad (Z 



To(z))(x,y)\ 



\z\"i S d g (x,y)~ 



n-l 

2 
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By Young's inequality, this provides the trivial bound 

n-3 111 

IKTHadO) -T (z))u\\ Loa ( M ) < |z| V|| u || Lp(M)) - < - + — 

and by interpolation, we obtain the following additional bounds 

n — 3 _ _n_ 

IKZkad - To)(z)u\\ Lq{M) < C\z\ * 2, ||n|| iP(M) 

.,1 n-11 - 1 1 1 

if - < : and - < - H . 

g n + 1 p' p 2 2n 

Since ?b(z) satisfies similar estimates with better decay in |z| when 1/p < 
1/2 + l/2n, the previous estimates remain true for the Hadamard para- 
metrix. These estimates with their dual provide the following additional 
boundedness properties of the Hadamard parametrix. 

Lemma 4.5. Let l<p<2<qbe such that 1/p—l/q < 2/n, the Hadamard 
parametrix satisfies the following bounds 



|rHad(^)IU(L9,LP) < C 




jf I < n=l J_ nnt j I < 1 _i_ J_ 
V g — n+1 p' a ' ia p ^ 2 + 2n 

y 9 — n-l p' a ' tCi 9^2 2n 



The range of exponents covered by this theorem is the union of the two 
violet trapezia in figure [6j 




Figure 6. Extended set of admissible exponents 
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In particular, we have 

||2Had(z)u| 



(4.16) 



< 



A 4 IMIl 2 (m) 



\Tuad(z)u\\ L 2 



(M\ S, \z\ 4 \\u\\ 2n 



These estimates will prove useful in the final argument proving the resolvent 
estimates when dealing with the remainders. 

4.3. Remainder estimates. We now deal with the remainder term S(z): 
the results are similar than for the Hadamard parametrix (in fact the range 
of admissible exponents is larger) but the proof is simpler due to the lack of 
singularities since the kernel is supported away from the diagonal in M x M. 



Lemma 4.6. The operator S(z) defined in (3.6) is a bounded operator 
(4.17) S(z) : L p (M) -> L q {M) 

when 1 < p < 2 < q with a norm bounded by 



\\S(z)\\ c{LP , Lq) <C< 



71-1 

\z\ 4 

n — 1 r 

\z\ 4 

n— 1 

\z\ 4 



2q 

V t> 
n 



if l -< 



re-1 1 



if ^Et A < - < , - ■ 

J n+1 p' — q — n— 1 p' ' 

i/ 1 > 0±1 J, 
J g — re— 1 p 



n+1 1 



Proof. The kernel of ^(z) is uniformly bounded by |z| _1//2 because of (3.9), 
hence £2(2) is a bounded operator from L P (M) to L q (M) with norm bounded 
by |z|~ 1//2 . For this part of the operator we actually get better bounds than 
expected since \z\ > 6. We now concentrate on the part Si(z) of the oper- 
ator: because of Lemma 3.2 and the support of grad p y), the kernel of 
this operator looks like 



n-l 

\z\ 4 e 



-^ d ^ y) b(x,y,z) 

where b is a smooth function. 

When argz G [—n/2,n/2], Young's inequality provides a better bound 

n + 1 , n f 1 _ 1 A 

\\Si(z)u\\ Lq(M) < \z\ 4 

When argz ^ [— 7r/2, 7r/2], the bound follows from Young's inequality, The- 
orems 2.2 and 2.4, and duality and interpolation. □ 



-c\z\ l / 2 \t II 

e 11 II«IIlp(a/)- 



4.4. Resolvent estimates. Before actually turning our attention to resol- 
vent estimates, let us first patch our local constructions into a global para- 
metrix. We cover the manifold M with finitely many open geodesic balls Uj 
of radius e contained in charts of M. We consider a partition of unity 



1 = ^2 X P sn PPXjCUj 
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subordinated to the covering {Uj)i<j<j ■, and define 

T Had( z ) to be the local 

Hadamard parametrix constructed in section [3] using Xj( x ) ® Xj(u) as the 
cutoff function. By a slight abuse of notation, we keep Tr^^z) to denote 
the global parametrix 

T Rad (z) = Y J T^ d (z) 
and similarly for the remainder 

J 

S(z) = J2s (j) (z). 

3=1 

Having constructed those operators, we indeed get a global right parametrix 
(4.18) (-A g + z)T Uad (z) = Id+S(z). 

As a finite sum of operators, THad(^) shares the same continuity properties 
as the local parametrices, and this is also true of S(z). 

We are now ready to prove the resolvent estimates. We use ideas coming 
from |14j and also [21 [3] to deal with the remainder terms. 

4.4.1. L? estimates. Let = Ao < Ai < • • ■ < Afc < . . . be the sequence of 
eigenvalues of the (positive) Laplacian — A g repeated with multiplicity and 
(V'fc)fcLo be t ne corresponding sequence of unit eigenfunctions. We denote 



f{k) = / fl/) k dVg 

Jm 

the Fourier coefficient of the function / corresponding to the eigenfunc- 
tion ipi~. 

Lemma 4.7. Let (M,g) be a compact Riemannian manifold (without bound- 
ary) of dimension n > 3, for all z £ C \ R_ we have 

IM| L 2 (M) < \z\~^{He^y l \\{-/\ g + z)u\\ L 2 {M) . 

Proof. We start with the factorization 

Afc + z = {\f\k +i\fz) (y/Xk~iy/z) 

which implies 

+ > (A fc + \z\)* Re y/z. 
If / = (— A g + z)u, using Fourier series, we get 

»*(*) = E^fTl2 ^ \znReVz\- 2 \\f\\h iM y 
This completes the proof of the lemma. □ 
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Remark 4.8. When Re \fz > 5, we deduce the following improved L 2 
bound 

( 4 -!9) \\Al 2 {M) ^ \z\~M\(- A g + z)u\\ L 2 {M) . 

which will prove to be of crucial importance to estimate the remainder terms 
in the proof of Theorem Unfortunately, such an improved bound only 
holds when the spectral parameter is outside of a parabola. Indeed suppose 
that the estimate ( 4,19[ ) holds then testing it with unit eigenfunctions ipk 
implies 

+ z\ > c\z\ 2 

and the parameter z has to remain outside of the envelope of the curves 
\Xk + z\ = c\z\ 2 , k £ N. In particular, if we choose z = —\k + iy then we 
have 

_ y*>c\\l + y% 

and \y\ > yJcXk- 

Besides, if one is interested in proving resolvent estimates when z lies in 
5,5, then difficulties arise when the parameter z lies on the lines Imz = ±<5. 
This means that the square root of z lies on the part of hyperbola 

r- x 1 6 , /2 
y/Z = X +1-X, A > \ - 
2 V o 

(see figure [7]) and this implies 

\z\~5 1 Re yfz\ 1 ~ 1. 
In other words, the I? estimate is weak 

\\u\\l.2(m) ^ \\{~\ + z )u\\l*{m) 

in terms of the decay with respect to the size of the spectral parameter 
Ul ~ A 2 . 



In In 2n 2n 

4.4.2. L n + 2 — L»- 2 estimate. We now prove the L n + 2 — L n ~ 2 estimate of 
Theoremll.il We start from 



u = T Ha d(^)(-A 9 + z)u + v 

where 

v = (ld-T Uad (z)(-A g + z))u. 
Note that by ( |4.18| ) we have 

(4.20) (-A g + z)v = -S{z)(-A g + z)u. 

We start our estimations 

IMI g»„ , < \\Tna.d(z)(-Ag + z)u\\ _2n +\\v 
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Figure 7. Values of the spectral parameters z G and 
and by the boundedness properties of the Hadamard parametrix given in 

Ag + Z)U\\ 



Theorem 4.1, we have 
llu" 



2n < 



2n + U 2n 



To estimate the second right hand-side term, we need the following lemma. 

Lemma 4.9. Let (M,g) be a compact Riemannian manifold (without bound- 
ary) of dimension n > 3, and let 5 be a positive number. For all z £ Hj we 
have 



(4-21) \\v\\ L ^ (M) <C\z\-4\\(-A g + z)v\\ L 2 {M) . 

Proof. We use the adjoint of Hadamard's parametrix as a left parametrix 

v = T HsLd (z)*(-A g + z)v - S(z)*v 
and get the following estimate by the boundedness properties of the Hadamard 



parametrix given in Theorem 4.4 — more precisely (|4.16|) — , and of the re- 
mainder given in Lemma 4.6 



< 



'2,n=2(Af) 



z)v\\ L 2 (M) + \z\±\\v\\ L 2 {M y 



The proof of (4.21) is completed once one has controlled the second right- 



hand side term using the bound obtained in Lemma 4.7 
IMIl2(m) £ \z\'^\\(-A g + z)v\\ L 2 (M) . 
This ends the proof of the lemma. 



□ 



Using the previous lemma and the identity (|4.20|), we get 

< 



\U\\ 2n 



Ag + Z ) U W L ^, M) + M ^ \\S{Z){-Ag + Z) U \\ L 2 (M) . 
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It is now a matter of using the boundedness properties of the remainder 



S(z) given in Lemma 4.6 to conclude that 



lull 2n < II ( — A„ + Z)u\\ 2n 



This completes the proof of Theorem 1.1 



5. Relation with LP Carleman estimates 

The aim of this section is to illustrate the relations between resolvent esti- 
mates and analogues of the LP Carleman estimates obtained in the Euclidean 
case by Jerison and Kenig |llj (for logarithmic weights) or by Kenig, Ruiz 
and Sogge (12] (for linear weights). The proof of L 2 Carleman estimates 
for limiting Carleman weights in [5] is based on integration by parts and 
can not be used in the LP setting. However in [3], we were able to prove 
such estimates following an idea of Jerison [10], (see also [THl Section 5.1]) 
using spectral cluster estimates of Sogge [E] . Here we present an alternative 
argument following Kenig, Ruiz and Sogge p2] and Shen [IT] and reducing 
the Carleman estimate to resolvent estimates. 

Limiting Carleman weights on an open Riemannian manifold (N, g) are 
smooth functions (p such that dip /OoniV and 

(5.1) {pj,P v }(x,Q = 0, (x,0€p- 1 (0)cT*iV 

if p<p denotes the following symbol 

Plp (x,o = \a\ 2 g -\dip\ 2 g +2i(d<p,o g . 

This a conformally invariant property. This notion was introduced in |13] 
which dealt with Calderon's inverse conductivity problem with partial data 
in dimension n > 3; this article was concerned with the isotropic setting and 
took advantage of replacing classical linear weights by logarithmic weights 
(which are examples of limiting Carleman weights in the Euclidean setting) . 
The motivation for introducing such weights in the context of general Rie- 
mannian manifolds lies in the application to inverse problems |13} [5j H] where 
solutions to the Schrodinger equation with opposite exponential behaviours 
are constructed on manifolds using complex geometrical optics. 



It was proved in OUS] that the condition (5.1) is in fact stringent. Sim- 
ply connected Riemannian manifolds which admit limiting Carleman weights 
are those which are conformal to one having a parallel field (i.e. simulta- 
neously a Killing and a gradient field). In particular, such manifolds are 
locally conformal to the product of an Euclidean interval and another open 
Riemannian manifold. The results in P H| are concerned with the sim- 
pler situation of manifolds conformal to global products of the Euclidean 
line with another Riemannian manifold. This is the situation in which we 
choose to set ourselves to underline the connection between resolvent and 
Carleman estimates. 
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The setting is the product N = RxMq of the Euclidean (one dimensional) 
line and of an (n — l)-dimensional compact Riemannian manifold (Mo, go) 
without boundary, endowed with the product metric 

g = dxj + g . 

We consider the Laplace-Beltrami operator on N 

P = A g = d 2 Xl +A go . 

A natural limiting Carleman weight on this product manifold is x\ and the 
corresponding conjugated operator reads 

e r^p e -r Xl = _ 2t Q xi +T 2 + Aso _ 

Theorem 5.1. Let N = R x Mo be the product of the Euclidean line and of 
an (n — 1) -dimensional compact manifold (Mo, go) endowed with the product 
metric g = dx 2 + go . For all compact intervals I C R, there exist two 
constants C, to > such that for all r G R with \t\ > To, we have 

(5.2) \\e TXl u\\ 2n < C1|e T:ri Pd| 2n 

when u G C X '(N) is supported in I x Mo. 

This Carleman estimate was obtained in [3], and in \17\ when Mo = T" _1 
is the n — 1 dimensional torus, by a different method, which uses as well 
the spectral cluster estimates of Sogge but on n — 1 dimensional manifolds. 
By translating and scaling the estimate in the first variable, it is always 
possible to assume I = [0, 2n] without loss of generality, which we do since 
this will lighten our notations. The proof follows [12] and is based on the 



construction of an inverse operator of the conjugated operator (5.3); to ease 
such an inversion, inspired by [8] (see also [El Section 3.2]), we further 
conjugate the operator by a harmless oscillating factor 

(5.3) _ e ™i-i*ip e -™*+!*i = (d xi + £) 2 + 2it(d xi + - r 2 - A go . 

We regard functions u G Cg°((0,27r) x Mo) as smooth periodic functions in 
the first variable. This corresponds to a change of perspective by replacing 
the manifold = RxMo by the product M = TxMo of the one dimensional 
torus T = R/27rZ and of the manifold (Mo, go)- We denote by Ao = < 
Ai < A2 < • . ■ the sequence of eigenvalues of — A go on Mo and (V>fc)fc>o the 
corresponding sequence of eigenfunctions forming an orthonormal basis of 
L 2 (M ) 

We denote by ix k : L 2 (Mo) — > L 2 (Mo) the projection on the linear space 
spanned by the eigenfunction ipk so that 



^2 TTfc = Id, X k TT k = -Ag Q . 

k=0 k=0 
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The eigenvalues of the Laplacian A g on T x Mo are — (j 2 + Afe) with j G Z, 
/c G N and the corresponding eigenfunctions ar^] 

We denote by ttj^ : L 2 (M) — > L 2 (M) the projection on the linear space 
spanned by the eigenfunction e^^ipk'- 

«27T 



Kj,kf(x) = ^(J e %jyi nkf(yi,x')dyi) e 

and define the spectral clusters as 

Xm= ^ Kjfa m G N. 

Note that these are projectors Xm = Xm- We end this paragraph by recalling 
the spectral cluster estimates of Sogge [HI [19] that we will need: 

\\Xmu\\ 2» < C(l + m)2||ti|| L 2 (M) 

(5.4) L "- 2(M) 

\\XmU\\ L 2 iM ) < CCl + mJaH^ . 

The two inequalities are dual to each other and can be found in [18, Corollary 
5.1.2]. 

We are now ready to reduce the proof of Carleman estimates to resolvent 
estimates. Recall that our goal is to prove 

(5.5) llull 2 „ <C||/|| 2 „ 
when u G Cq°(( — 7r, 7r) x ilfo) and 

(5.6) (d xi + ^ u + 2ir (yD Xl + n - A 9o u - r 2 n = /. 

The equation (5.6) is actually easy to solve: writing / = J2jk n j,kf an d 
similarly for u, the equation formally becomes 

j + ^ + 2iT (j + - T 2 + Afc^j TTj- fcU = 7Tj- jfc / 

for j G Z and for k G N. The symbol on the left is always nonzero provided 
that \t\ > 1. Thus, the inverse operator may be written as 

OO OO j, 

q f = V" y 7 [Ml 

T ^octoO + h) 2 + Hj + h)T + Xk-T 2 ' 



Here we are using in an essential manner the fact that the metric g is a product, i.e. 
the fact that x\ is a limiting Carleman weight. 
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We now use Littlewood-Paley theory to localize in frequency with respect 
to the Euclidean variable x\\ we have 

oo 

U = 

with 



u v = E ^( r e~ ijyi u{ yi ,x')dy^\ji x \ u>0 
and similarly for /. It suffices to prove 

(5.7) ||liJ| 2n < CII fjl 2n 

since then by summing the former inequalities and using Littlewood-Paley's 
theory 

oo oo 



L"- 2 *■ — ' L"~ 2 * — ' L™- 2 L"+2 

i/=0 ^=0 



and this is the inequality (|5.5|) we were looking for. 



We now concentrate on the microlocalized estimate (5.7). Because the 
conjugated operator (5.3) and the localization in frequency l[2"-i 2")(Aci) 
commut^J we have 

u v = G T fv 

We denote 

2 \ -1 



R(z) = ( ( D X1 + - ) - A 90 + z 



r 

i\ 2 



the resolvent associated with the elliptic operator (D Xl + |) — A 9o . In the 
case v = 0, we have 

(5.8) n = G T / = R(-t 2 + ir)/ 

so ( |5.7[ ) is a resolvent estimate with 2 = — r 2 + it. When v > 1, we want to 
evaluate on microlocalized functions the error made by freezing the operator 
2tD Xi into 2 u t in (5.3). This amounts to replacing G T by the resolvent 
R(—t 2 + i{2 v + 1)t). More precisely, we want to estimate the term 

(5.9) (R{-T 2 + i{2 v + \)r)-G T )f v . 
This can be explicitly computed 

oo oo 

( R{ _ T 2 + i{2 » + i)T)-G T )u= E EW^/- 

j=— oo fc=0 



4 Here again we are using the fact that the weight xi is a limiting Carleman weight, i.e. 
that the metric go does not depend on xi. 
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where 
a U jk( T ) 



ir(2--2j)l [2 ,-i )2 „ ) (i) 



j + \f + 2ir{j + i) - r 2 + A*) ( (j + \f + i(2- + 1)t - t 2 + X 



Using the spectral cluster estimates (5.4) and the fact that spectral clus- 
ters are projectors, we have the following string of estimates 

\\(R{-t 2 + i{2 v + l)T)-G T )f u \ 



2n 



< C 1 + "0* \\Xm{R(-T 2 + i(2" + l)r) - R r )f u \ 



m=Q 

oo 



L 2 (M) 



<^(l+m)2 SUp \a U jk (T)\ X \\x m fu\\L^(M) 

m=0 m<\J j 2 +\k<m+l 



and using again (5.4), we get 



||(i?(-r 2 + i (2- + l)r)-G r )M| L ^ (M) 

<(X)(l + m) sup |< fc (r)|) x ||/, 

^m=0 m<y'i 2 +A & <m+l ' 

We now prove that the above series converges and is uniformly bounded 
with respect to r and v; we have 

2!/ | T l 

sup \ a ^)\ < { 2 _ 2)2 +1 2 

as well as 



^ 2"|r|(l+m) < /- 00 2 1r|t 

Cm2 _ t-2^2 i 4i/+1t-2 ~ / f+2 _ ~2\2 i Aiy+l^ al 



( (m 2 - r 2 ) 2 + 4^+V 2 ~ J Q (t 2 - r 2 ) 2 + 4^+V 2 

and if we perform the change of variables s = A~ v ~ l T~ 2 (t 2 — r 2 ) in the 
right-hand side integral, we obtain the bound 

^ 2 v \T\{l + m) < f 00 ds 

( m 2 _ T 2) 2 + 4 u+i T 2 ~ y_ oo 

Summing up our computations, we have the error estimate 

\\{R(-T 2 + i(2 U + \)t) -G T )fJ\ 2„ < HO 2n 

this means that in order to obtain ( |5.7[ ), it is enough to prove the resolvent 
estimate 

||i2(-T 2 + i(2 w + lYr)/J| 2 „ <\\U\ 2™ 
with a constant which is uniform in r and z/. 
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The conclusion of these computations is that Carleman estimates reduce 
to resolvent estimates of the form 

2 



(5.10) 
with 



U 2r, 



< 



Di + 



1 



Ago + z ) u 



2n 

LW+2(M) 



-r 2 + IQT (g > 1). 



After changing u into e z Xl u in (5.10), the resolvent estimate is a conse- 
quence of Theorem |1.1| since 



Re \fz 

provided |t| is large. 
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